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Abstract 

We analyze a class of quantum operations based on a geometrical representation of 
d— level quantum system (or qudit for short). A sufficient and necessary condition 
of complete positivity, expressed in terms of the quantum Fourier transform, is 
found for this class of operations. A more general class of operations on qudits is 
also considered and its completely positive condition is reduced to the well-known 
semi-definite programming problem. 
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1 Introduction 



Almost all quantum systems in the real world are open in the sense that they 
suffer from unwanted interactions with the outside world. The dynamics of 
an open quantum system is usually much more complicated than that of a 
closed one. One of the most important mathematical tools describing the dy- 
namics of an open quantum system is quantum operation, which has been 
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systematically studied after Kraus's seminal work [1]. Comparing with other 
mathematical formalisms coping with open quantum systems, such as mas- 
ter equation, Langevin equation and stochastic differential equation, quantum 
operations are especially appropriate for depicting discrete state change in 
discrete time. Recently, the rapid development of quantum information pro- 
cessing technology revives a wide interest on quantum operations due to the 
fact that quantum information processing systems, for example, quantum com- 
puters, suffer from outside noises inevitablly, and that mainly a discrete-time 
evolution is concerned in these systems. Indeed, quantum operation plays an 
important role in many active fields such as quantum computation, quan- 
tum information, quantum error-correcting codes and quantum fault-tolerant 
computation [2,3,4] (for an excellent exposition, see [5] Chapter 8). 

There are several equivalent ways to introduce the notion of quantum opera- 
tion; one of which is given in terms of completely positive (often abbreviated 
as CP) mappings. Let Ti. denote the Hilbert space of the principal quantum 
system. A linear mapping k on ?i is positive if it always sends a positive op- 
erator to another positive one; and k is completely positive if, furthermore, 
the mapping I ® k is also positive where / is the identity mapping acting on 
an arbitrary ancillary system. Complete positivity is a natural requirement of 
the real physical world without which the state of the composite system may 
be invalid after operations on its subsystem. Thus, a quantum operation is 
defined as a completely positive and trace-preserving linear mapping on the 
state space Ti. Some representations of the CP mapping have been presented 
in [1], also in[6], and an important one is the operator sum representation: 

N 

m-Y.^kpAi (1) 

where A^. is the conjugate transpose of A^. To guarantee the trace-preserving 
property, the following completeness condition 

N 

Y^AAk^I (2) 

is also required. For some other descriptions of quantum operations, we refer 
to [5]. As trace-preserving quantum operation is a fundamental mathemati- 
cal description of quantum channel, a thorough study on quantum operation 
will help us to understand the limits or capabilities of quantum information 
processing. 

The main goal of this paper is to find the necessary and sufficient conditions 
for certain kinds of linear mappings to be quantum operations. To motivate 
our problem, let us consider the trace-preserving quantum operation k on a 
two- level quantum system (or qubit for short). We employ the very useful 
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geometric tool of Bloch sphere to represent a qubit. Then every qubit can be 
depicted as a vector in it. Recall that the general form of the affine mapping 
on the Bloch sphere is 

r = (r^;, Ty, Tz)^!^ ^ M{r^, ry, + (C^;, Cy, Cz), (3) 

where r and are respectively the original and the image vectors under the 
mapping, M is a 3 x 3 real matrix and c = (c^., Cy, Cz) the displacement vector. 
It is easy to see that the effect of a quantum operation on qubits is just an 
affine mapping on the corresponding vectors. Conversely, a question naturally 
arises: whether or not every affine mapping on the Bloch sphere has a cor- 
responding quantum operation? The answer to this question is unfortunately 
no. In fact, there are very simple affine maps that cannot be images of quan- 
tum operations [7,8]. A sufficient and necessary condition of the special case 
without any displacement has been found by several research groups [7,9,10] 
and some special cases of {c^, Cy, c^) are also considered. In [11], M. B. Ruskai, 
S. Szarek and E. Werner have completely solved the general problem for the 
qubit case. More precisely, given the transition matrix M and displacement c, 
the sufficient and necessary condition of when such a mapping is completely 
positive was presented there. 

A qubit is the simplest quantum system which has been studied thoroughly. 
However, quantum systems of higher dimension also concern us in quantum 
information processing and are relatively less studied. This observation leads 
us to consider higher dimensional generalization of the above question. To the 
authors' best knowledge, the higher dimensional version of this question is 
still open. In the present paper, we will carefully examine the affine mappings 
in higher dimensional quantum systems, and give a partial answer to this 
question. This will certainly give some new insights to the field of quantum 
information processing and bring a useful mathematical tool for the solutions 
of some problems in this field. 

The paper is organized as follows. In Section 2, we extend the vector repre- 
sentation from qubit to qudit case. This enables us to give a (i-dimensional 
generalization of affine mappings. Some simple properties of such geometric 
representation are presented. In Section 3, we focus on a special class of higher 
dimensional affine mappings where the transition matrix is diagonalized and 
the displacement vector vanishes. A sufficient and necessary condition of when 
an mapping in this class is a quantum operation is found. With this condition, 
one can easily check whether a given mapping in the special case is an image 
of a quantum operation or not. Some well known results, for example, the uni- 
versal NOT gate with optimal fidelity in the high dimensional state space, can 
be recovered. A more general case involving multiple qudits is also discussed. 
The most interesting thing is that this condition is deeply related to the quan- 
tum Fourier transform (QFT for short), and such a surprising connection is 
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not easy to observe if we only deal with the case of qubits. In Section 4, we 
further consider higher dimensional affine mappings with diagonal transition 
matrix in which the displacement vector presents. They are of the most com- 
mon affine mappings in practical use. The problem of determining when such 
an affine mapping is a quantum operation is then reduced to the semi-positive 
definite programming problem. We argue that this problem has no simple so- 
lution in the general case. However, we are still able to give a simple sufficient 
condition of it. Finally, we draw a brief conclusion in Section 5. 



2 A geometric representation of qudits 



It is well known that for qubit system, a wonderful geometric representation 
using Bloch sphere exists. To be specific, the density matrix p of a qubit can 
be written as 

I -\- r ■ a 



where r = (r^;, r^, r^) is a 3-dimensional real vector with Euclidean norm | |r| | < 
1. a = (X, y, Z) is a formal vector of Pauli matrices and f-a = rxX+ryY+rzZ 
where X, Y, Z are Pauli matrices. 
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It is notable that ||r|| = 1 if and only if p is a pure state, and p is a mixed 
state whenever ||r|| < 1. Given a vector r, we can completely determine the 
corresponding quantum state p and vise versa. Furthermore, the dynamics of 
p can be visualized through the affine operation on f. The key of a qubit pos- 
sessing the Bloch sphere vector representation is there exist a set of matrices 
{X, Y, Z} which, together with the identity matrix /, forms a basis of the 2x2 
matrix space. 

To extend such a representation to a higher dimensional state space, we choose 
a set of generalized Pauli matrices [2] as the basis. Suppose we have an or- 
thonormal basis {\k)} of a qudit space 7i, if we define X a.s X\k) = \k + 1), 
where the addition is modulo d, and Z as Z\k) — u'^lk), where oj — exp{—^) 
is the dth unity root, then the whole set of the general Pauli matrices in qudit 
system can be expressed as {(Tp^q = X^Z'^ : p,q = 0, ■ ■ ■ ,d—l}. For a detailed 
presentation of the properties the generalized Pauli matrices enjoy, we refer 
to [2] . Here we only list some of them related to our results (for simplicity, we 
define Xi ^ X\ Zi ^ Z'): 
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(1) (Multiplication relations) 



(2) (Hermitian relations) 



(3) (Commutation relations) 



= ^'^<^{d-i),{d-3)- 



(4) (Orthogonal relations) Suppose the inner product ol dxd linear space is 
defined as (A, B) = tr(AtS), then 



(<7jj, o'i'j') — d6i^i>6jjr, 



where S is the Kronecker function. 



1, ^/ i = 3 



(5) The Pauli matrices {ai j : i, j = 0, 1, ■ ■ ■ , rf — 1} form an orthogonal basis 
oi d X d linear space Q. That is, for any ^4 e Q, 



d 

The next proposition shows that on a higher dimensional space the matrix 
composed of all eigenvectors of X corresponds exactly to the Fourier transform 
F. 

Proposition 1 Suppose F is the quantum Fourier transform on d- dimensional 
Hilbert state space H, that is 

F\j) = ^j:e'-^^'/'\k), 
yd k=o 

then Xk — F^ZkF, where {\j) : j — 0,1,2, ■ ■ ■ ,d— 1} is an orthonormal basis 
ofH. 

Proof. Simple, details omitted. □ 

Using this set of generalized Pauli matrices in high dimensional space, we can 
decompose a density matrix p of a qudit in the principal Hilbert space Tl as 
follows: 

— * — * 

r • a 
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where f = {tr {a qq p), tr {a^i p), ■■ ■ , tr((T^_]^ is a d^— dimensional vector 

and a is the formal vector such that a — (ctq^o, cro,i) ' ' ' ■> (^d-i,d-i), the formal 
inner product is then defined as 

d-l 
p,q=0 

The representation here is compatible with the Bloch vector notation since we 
can simply add / and 1 to the Pauli matrices a and f respectively. Moreover, 
some properties of Bloch vector notation preserve in the new representation. 
For example, = Vd when p is a pure state; < y/d when p is a mixed 
state, while ||r|| = 1 if and only if p is the maximally mixed state I/d. It is 
also useful to point out that to make p a legal quantum state, we always have 
ro,Q = tr{aloP) = tr{p) = 1. 

Suppose quantum operation £ onTi maps p to p' with the representation 

a simple computation tells us there exists a x dimensional matrix M 
such that 

r' = Mr. 

In other words, every quantum operation leads to an affine mapping with the 
form 

f^Mf. (5) 

Then the question raised in the introduction can be more explicitly stated as 
follows: is every affine mapping with the form (5) a quantum operation on state 
space 7i? In general, the answer to this question is certainly no. Our purpose 
of this paper is then to characterize the affine mapping that corresponds to 
some quantum operation. 

Before we go deep into this question, let us first examine some properties that 
transition matrix M and vector f must satisfy. We say a quantum operation 
S is unital if S{I) = I. Otherwise we call it is non-unital. Suppose transition 
matrix M corresponds to a trace-preserving quantum operation £, then M 
has the block form 

A 

M = 

\c M'^ 

where c is a c?^ — 1 dimensional vector (we often call it as displacement vector), 
M' is a (o?^ — 1) X [d"^ — 1) matrix. Furthermore, if S is unital, we have c = 
(vector). For p is a positive operator with trace one, we get f- a = (f ■ a)^ and 
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Too = 1, or equivalently, 

^P,q = u}'^'^r-p-q,roo = 1, 

where all the arithmetical operations of the indices are modulo d, so V-p-q ~ 
^d-p,d-q- The above equations can be considered as a general constraints on 
the vector r. With this condition, we can use — 1 real numbers to represent 
a qudit state p. 

Denote the set of vectors satisfying the above equations as V^2. A simple 
calculation shows that \/^2 is close under addition and multiplication by a real 
number; that is, rl + r2 G V^a and tfl G Vcp provided that ri,r2 G Vcp and 
f is a real number. Notice that when i is a complex number, tf{ G Vip. does 
not necessarily hold. Suppose A is a x matrix which is a linear mapping 
on V(pi then Af G V^p for all r G V^2 if and only if {Ar ■ = Ar ■ a for any 
r G V^2 . Let us further consider a special but very important case where A is 
a diagonal matrix 

A = diag{\o, Ao,i, • • • , \d-i,d-i), 

or simply we write A = diag{X). Then the condition turns out to be A_p^_g = 
A*^^. So if a quantum operation maps p with the Bloch vector r to S{p) with 
the Bloch vector Mr, where 

M = diag{Xofl, Ao,i, • • • , Xd-i,d-i) = diagiX) 

is diagonal, then the following conditions must be satisfied: 

fp,q — <^ T_p_q, '^-p-q — \^q- 

For the case of qubits, i.e., d = 2, we have the pauli matrices 

c^o,o = I, cro,i = cTi^o = X, ai^i = XZ 
and cu — —1, then the above conditions about f may be rewritten as 

ro,o = 1, ro,i = r*_i, r^o = rI,o,^i,i = -rli, 

that is, ro,o = 1, '"o,!) ''1,0 ^-re real numbers and ri^i is pure imaginary number. 
The conditions about M can be rewritten as 

Ao,o = 1, Ao,i = Ag 1, Ai,o = Ai Q, Ai,i = Xli, 

all the four numbers ^0,0, ''"o,i) ''"1,0 and ri^i are real. Moreover, by introducing 
Y = iXZ, we can also make ri^i a real number. Then the qubit state p and 
real vector r = (ro,o, to,i, ^"i,o, ''"1,1) are one to one corresponding to each other 
when |rj < 2, and the scalars of the four axes are all real numbers. 
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3 Afflne mappings with diagonal transition matrix 

In this section, we consider the special affine mappings with diagonal transition 
matrices M and no displacement vectors. The effect of these mappings is just 
a scalar multiplication by Xp^g in the direction ap^g, or 

r I— >■ Mr, (6) 

— * 

where M — diag{\) is a diagonal matrix, we can also write 

^ {'^P,q) = K,q(^P,q^ P, ? = 0, • • • , d - 1. 

It is well known that when qubits are considered, with a suitably chosen 
axis, many interesting quantum operations have this simple form in the Bloch 
sphere. The typical examples are bit flip channel and phase damping [5] . 

This section is devoted to find the necessary and sufficient condition of such 
affine mapping to be a quantum operation's image. As mentioned before, a 
key feature of a quantum operation is its complete positivity. To determine the 
complete positivity of we need the following theorem attributed to Kraus 

[1]: 

Theorem 1 Suppose S is a linear operation on Hilbert space H, IZ is an 
ancillary Hilbert space with the same dimensionality with H. Let {\k)} be an 
orthonormal basis for H and TZ. Then £ is CP mapping if and only if the 
operator 

{I®£){\a){a\) 

is positive, where 

\a)^Y.m)l^d 

k 

is the maximal entangled state of the composite system. 

Using Kraus's theorem, to decide whether a given linear mapping is a CP 
one, wc need only examine the positivity of eigenvalues of the operator (7 (g) 

E){\a){a\). 

We now are able to present one of the main results of this paper, namely, a CP 
condition for an affine mapping E with a diagonal transition matrix and no 
displacement vector. To one's surprise, this condition is essentially connected 
to the quantum Fourier transform F, which, as we all known, plays a crucial 
role in quantum computation (for example, quantum Fourier transform is a key 
step for the Shor's Factoring algorithm and Discrete Logarithm Problem[4]). 
To be specific, we have the following theorem. 
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Theorem 2 The linear mapping defined in (6) is a CP one if and only if 

{F®F^)\ > 0. 
where X = (Ao,o, Ao,i, • • • , Xd-i,d-i)- 

Proof. Using the Kraus's theorem, we need to calculate {I ^ £){\a){a\). So 
first we expand \k){l\ as the linear combinations of {cij}- By the orthogonality 
of (7ij, we get 

m,n 

where am,n = tr{al^Jk){l\)/d = uj-~'^^8m+i,k/ d. So 
{I®£){\a){a\) 



^Y.\k)m\k){l\)/d 

k,l [I) 

= ^~^'^^"'^^'^m,g^Crm,nXm,n/d^ 

l,m,n,q 

m,n 

The crucial point of the problem is that the set 

form an Abelian group of order d^, which means that this set is closed under the 
operation of matrix multiplication, matrix inversion, and its elements commute 
with each other. We prove commutation relation as follows: 



(8) 

Thus, we can simultaneously diagonalize the set of matrices {(Tp-q ® CTp,?} as 

(7m -n ® (^m,n ^UD{m,n)U\ 
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where D{m, n) is a diagonal matrix, and U a unitary matrix. It must be noted 
that the matrix U can be chosen as the same for all the matrices in the set. 
Then 

^ m,n 

The diagonal entries give the desired eigenvalues of (/ <Si S){\o£){a\). Define 

one can easily check 
So by (7) we have 

(I®£)(\a){a\)\^s,t)^l^s,t\^s,t), 

where 

^l/dJ2{F)s,m{F\^\^,n (9) 

m,n 

= lM(F®Ft)A],,t, 
or in matrix representation, 

/2= l(F®Ft)A. 

For (/ (g) S){\a){a\) is a Hermitian matrix, it is positive if and only if all its 
eigenvalues are non-negative. This leads to the relation 

(F(8)F^)A > 0. 

With this we complete the proof of the theorem. □ 

There are many ways to gain the result in the theorem and some of them are 
simpler. But the method given here is a coherent one in the sense that when 
the displacement vector does not vanish it still works. 

Now we consider further the case where N qudits are involved. Suppose 

the generalized Pauli matrices on rf-dimensional Hilbert space Ti. are g : 
p,q = 0, ■ ■ ■ ,d — 1}. Use this set Pauli matrices, we can obtain an orthogo- 
nal basis of operators acting on 7Y®^. For the sake of convenience, we denote 
i — (^l,^2, ■■■j^at) and j — (ji, j2, ■■■,jN), then we define the orthogonal basis 
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as TTij = (Til J, (g) (Ti2j2 ■ ■ ■ , ,where • • • ,iN,jN are range from 

to d — 1. Then the following theorem is a straightforward generalization of 
Theorem 2. 

Theorem 3 Let H be a d- dimensional Hilbert space, and let S be a linear 
mapping on Hilbert space such that 

Then S is CP if and only if 

(F®Ft)®^A> 0. 

Proof. Same as the previous theorem. We omit the details here. 

Wc should point out here that the quantum Fourier transform F in the above 
theorem is defined on 7i. There is another Fourier transformation F' which is 
defined on A^d-dimensional state space 7Y®^. F' is not equal to F®^ in general. 
This enables us to obtain two different conditions. However, this two conditions 
are not related since they describe different types of affine mappings which are 
directly related to the choice of basis. One basis is general Pauli matrices and 
the other is tensor products of general Pauli matrices. 

As a conclusion of this section, we show some simple applications of the above 
two theorems. First, let us come back to the case of qubits. Notice that in this 
case, quantum Fourier transform F is just the Hadamard gate 

H={X + Z)/V2 = ^ 

For = H, the CP condition given in Theorem 3.3 can be rewritten in a 
more brief form: 

H^'^^'X > 0. 

When only one qubit involves, the condition deduces to if®^A, or more explic- 
itly (noticing that Ao,o = 1) 

1 + Ao,i + Ai,o + Ai,i > 0, 

1 ~ Ao,i + Ai^o ~ Ai,i > 0, 

which is similar as the equation (12) in [11]. 

Since {F ® F'I")®^A has an entry J2^i,j/d'^ , a necessary condition for £ to be 
a CP mapping is that J2 \,j > 0. An interesting special case here is when all 
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entries in A have some same real value p except Ao,o = 1- This is in fact a gen- 
eralized depolarizing channel. In this case, we have two different eigenvalues: 
1 + {(f — l)p and 1 — p. The CP condition then becomes 

-l/{d^ - 1) <P < 1. 

This result is well known and here it comes out to be a direct corollary of our 
above theorem. 

What we would like to point out here is that when c is negative , — l/((i^ — 1) 
is the limit value we can achieve in this operation. It is also easy to see that 
this value decreases with d increasing. This means that when dimensional- 
ity increases, the state space becomes more complicate, and the constraints 
between the different axes are more rigor. 

Let us consider further the case where p — — — 1). The pure state p — 
will be mapped to £{p) = {dl — p)/{d'^ — 1). If we treat the operation 
as an unprecise universal NOT operation, then the fidelity will be 

m£{p)m = (10) 

which has been proved optimal in the case of c? = 2, see [12,13]. 

Recall from Section 2 that we have the constraints Xp^q — X*_p _g. So these 
operations are not completely independent in different axes, and there are 
conjugate axes such as CTp^g and (y_p _q. Thus, the scalars of a CP operation 
in this two directions must be conjugated with each other. In the qubit case, 
since are all self-conjugated, the entries of A are all real numbers, and the 
operation is independent for different axes. 



4 A more general class of afHne mappings 



What concerns us in the previous section are affine mappings without dis- 
placement vectors. In this section, we consider a more general affine mapping 
of the form: 



where 



Mr, 



M = diag{\) + 



fill 
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We denote the first column of 



'^0 0' 



as c. Intuitively, we call c as displacement 

c oy 

vector. Such kind of affine mappings correspond the most useful non-unital 
quantum operations in practice. The problem in which we are interested re- 
mains when such an affine mapping is an image of a quantum operation. The 
following theorem is the another main results in our paper which establishes a 
connection between this problem and the semi-positive definite programming. 

Theorem 4 The linear mapping defined in (11) is a CP one if and only if 

I ®C- a + ^ \p,q(Jp-q ® Up^q > 0. 
P,<1 



Proof. Prom 
we get 
and 



^((7o,o) = Ao,o(7o,o + c- a, 
^ {crp,q) = \,qcrp,q, {p, q) ^ (0, 0) 



1 

{I ®S){\a){a\) = —[I®c-a + ^Xp,qap,^q®ap,q\. 

" p,q 

Then by the Kraus's Theorem, the claim in the theorem holds, and we com- 
plete the proof. □ 

Since the set of matrices {(Tp _q(8>cTp^g} can be simultaneously diagonalized, we 
can rewrite the condition in the above theorem as 

U\l0c-a)U + D> 0, 

where D = diag{{F O F^)X)/d. 

There arc many basic problems in quantum information fields which can be 
reduced to the most general version of the semi-positive definite programming 
problem; for examples, quantum state discrimination [14,15], state estimation 
[16], and quantum pattern recognition [17]. The above theorem provides us 
with an additional example of such problems. Indeed, given c. To keep £^ as a 
CP mapping, we need to determine the range of the diagonal matrix D. Then 
through D we can immediately get the range of A, for 

X^d(F^ ®F)D, 

— * 

where D denotes the vector formed by the diagonal entries of D. Purthermore, 
the problem of determining the range of a diagonal matrix D to keep A+D > 
for a Hermitian matrix A is exactly a semi-positive definite programming 
problem. 
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As is well known, there has not yet a general analytic sohition for the problem 
of semi-positive definite programming. We then cannot get a sufficient and 
necessary condition in an explicit form under which the linear mappings (11) 
are completely positive for the most general case. Here we propose instead a 
sufficient condition: 

Corollary 1 A sufficient condition for the linear mapping defined in (11) to 
he a CP one is 

M\< (^min, (12) 

where 

Afmin = mins,t[(-^ <S) F^)A]s,t/d, 

and 

is the norm of the displacement vector. 

The intuitive meaning of the conclusion is that the affine mapping (11) will 
be an image of a quantum operation if the norm of the displacement vector c 
is small enough. Of course, here c should satisfy the equation c - a — {c- a)^. 

Proof. Let c - a — VDcV^ be the spectral decomposition of c • ct, where Dc is 
a diagonal matrix and V is some unitary matrix. Since tr{c- a)^ — ||c|p, we 
have Ac < ||c||, whereAc is the maximal absolute value of eigenvalues oi c- a. 
So when (12) holds, /Xmin/ + Dc>0. Finally, we have 

I ® C - a + UDW > llrauJ + I <^ c ■ a 

^I®V{ii^iJ + D,)V^ (13) 
>0, 

where V is unitary. This completes the proof. □ 

The term I ® c • a -\- UDW is quite complicated, so one may hope to derive 
a simple expression of c, as what was done in the last section, which can be 
used to characterize the positivity of the original term I <S) c ■ a + UDU^. 
Unfortunately, this is impossible in general. The reason is that UDW and 
I <S> c ■ a do not commute with each other, but they are orthogonal to each 
other; in other words, we have 

tr{UDU\l ®c-a)) = 

This implies that we cannot diagonafize these two terms simultaneously. Oth- 
erwise, one of them will be led into a trivial case. Suppose that 

UDU^ = /, Ap,5 = 0, {p, q) ^ (0, 0), Ao,o = 1- 
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To keep 

/ (g) (/ + c • a) > 0, 

it only needs 

l + c-a>0. 

In other words, in this special case, the only condition that guarantees the left- 
hand side of the above inequality is positive is then not to move c out the unit 
sphere (of course, we need to require that c is a valid geometry representation 
of a quantum state). 



5 Conclusion 



In this paper, we use the set of generahzed Pauli matrices as a tool to present 
a new geometric representation of a quantum state of qudits. This enables us 
to introduce affine mappings into the space of qudits. Then we find a sufficient 
and necessary condition for a special class of affine mappings. This condition 
is also generalized to the case of N qudits. Using this result, we can easily 
decide whether a given mapping in the class we specified is a CP mapping or 
not. As an apphcation, we derive the optimal fidelity of the universal NOT 
on d-dimensional state space: 



Furthermore, we can easily deduce that when d = 2, or in the qubit case, this 
fidelity is optimal with the value 

The most interesting thing that we discover here is that the CP condition 
for this kind of affine mappings is closely related to the quantum Fourier 
transform. 

We also consider some more general affine mappings with diagonal transition 
matrices and displacement vectors. We point out that the problem of deciding 
the complete positivity of these mapping is equivalent to the semi-positive 
definite programming problem. So, it is impossible to derive a general symbolic 
solution of this problem. However, we still get a sufficient condition which says 
that when the norm of the displacement vector is small enough, operation 
under consideration is CP. 

Acknowledgement: We thank anonymous referee for helpful comments and 
informing us of a valuable reference. 
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